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, We propose a new kind of quantum phase transition in phase separated mixtures of Bose-Einstein 

condensates. In this transition, the distribution of the two components changes from a symmetric 
£\j | to an asymmetric shape. We discuss the nature of the phase transition, the role of interface tension 

and the phase diagram. The symmetric to asymmetric transition is the simplest quantum phase 
transition that one can imagine. Careful study of this problem should provide us new insight into 
this burgeoning field of discovery. 
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There is much recent interest in quantum phase transitions. Examples of these include the Wigner electron solid 
' melting transition, the Mott-Hubbard metal-insulator transition, and different magnetic transitions. In physical 
phenomena involving Bose-Einstein condensates (BECs) quantum mechanics play a crucial role. In this paper we 
propose a new kind of quantum phase transitions in phase-separated mixtures of BEC condensates. In this transition, 
the distribution of the two components changes from a symmetric to an asymmetric shape. To explore this transition, 
we first investigate the stability of the symmetric phase by studying its normal modes. We find interface modes that 
become soft. When the lowest frequency becomes zero, the instability sets in and this determines the stability limit 
of the symmetric phase. We determine the actual phase boundary by comparing the energy between the symmetric 
phase and the asymmetric phase and find that this actual phase boundary and the instability boundary is not the 
same. This suggests that the transition is first order. The system may be a good laboratory to study issues of quantum 
mctastability and tunnelling. The symmetric to asymmetric transition is the simplest quantum phase transition that 
one can imagine. Careful study of this problem should provide us new insight into this burgeoning field of discovery. 
We now describe the results in detail. 
q ' Mixtures of trapped Bose-Einstein condensates have recently received considerable theoretical [Q|2j|3|,[]]||[3]j7| and 
experimental interest. Experimentally at low fields, the spin exchange process can occur in an optically trapped 

condensate, leading to spin domains [[To) with metastable behavior [ ^HJl^l . Binary condensates in two hyperfine levels 
of 87 Rb have been created and studied |§Jl^], most notably realizing a system of interpenetrating Bose fluids ||, 
measurements of phase dispersion JlJ] , and a vortex state in a dilute-gas BEC [|l5) . 

The equilibrium density distributions of segregated mixtures in the absence of gravity have been studied numerically 
for different system parameters. Two types of configurations have been discussed: a symmetric HJH]]^] configuration, 
for which one component is inside the other one, and an asymmetric one in which the two components occupy the 
^vq ' left and the right hand side of a sphere . 

In general, the asymmetric phase possesses a lower interface energy. On the other hand, since the degree of 
self-repulsion may differ between the two species, the less self-repulsive component will prefer to remain where the 
density is higher, while the other component moves to the low density regions outside. This favors the symmetric 
phase. Depending on the system parameters, one of these two energetic considerations will win out. These system 
parameters can be adjusted by changing the trapping frequencies, the relative particle numbers of the two species, 
and the interaction between the particles with Feshbach resonances. We first address the stability of the symmetric 
O '. phase. 

O ■ We consider the two-component BEC in a spherically symmetric trap. The dynamics of this system is described 
^ , by the time dependent Gross-Pitaevskii equations: 

ih^ = -^Ay 1 + V ti y 1 + —(a 11 \* 1 \ 2 + a 12 \* 2 \ 2 )* u (1) 

■ at 2m m 

_ _ i 

iU ^T = "|- A *2 + ^tr*2 + ^-(a 22 |* 2 | 2 + ai 2 |*i| 2 )*2, (2) 
at 2m m 

where ^1,2 are the condensate wave functions, Vtr = muj^r 2 /2 is the trapping potential, u>q is the trapping frequency, 
r is the radial spherical coordinate, ay > are s— wave scattering lengths. 

We shall assume in this paper that the condition a\ 2 — 0-11022 > is satisfied and, therefore, the condensates 
are phase-segregated. We study the condensates in the TF limit. In this regime, the phase-segregated condensates 
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overlap over the length scale A = £/ \fa^]^Ja^a%i — T, where £ is the healing length |t]] . For parameters of the 
JILA experiments on phase-segregated states A « 47£. If the penetration depth A « ii, where R is the size of 
the system, the condensates can be approximately treated as nonoverlapping, which we assume to be the case. The 
effect of overlapping condensates results in a finite surface tension and can be included via boundary conditions at 
the interphase. If the condensates do not overlap one can neglect the last terms in Eqs. ([!]), (Q). As a result, the 
dynamical equations for \Pi and ^ 2 decouple. However, the two condensates are coupled by the boundary conditions 
at the interface which require continuity of the pressure and the normal velocity. 

The symmetric phase consists of a central core of the first component and an outer shell of the second species (we 
assume an < 022). The stationary density distribution m = \^i\ 2 of the two components is given by 

Gum =//i(l- r 2 /Rl), 0<r<R„ (3) 



G 22 n 2 = /x 2 (l - r 2 /R 2 2 ), R* < r < R 2 , (4) 

where Gu = Anf^au/m, Ri = y // 2fii/muj 2 . The normalization condition J riidV = A;, where are the numbers 
of condensate particles, determines the chemical potentials \Xi. The position of the phase boundary i?» is given by 
the condition that the pressures exerted by both condensates are equal R* = R 2 y/(1 — n\) /(l — k)A JTsJ , where we 
introduced the dimensionless parameters k = \/an/a 22 , X = fj, 2 / Hi- The symmetric configuration is favorable when k 
differs from unity, with the less repulsive component in the middle ( K < 1). At the interface between two components 
n 2 /ni = k < 1. 

One can rewrite the time dependent Gross-Pitaevskii equations in a hydrodynamic form which shows an analogy 
between our problem and motion of two immiscible fluids. In the strong phase-segregated regime the dynamics of 
each components is described by the following Josephson hydrodynamic equations |l7| 

^+V.(nV) = 0, (5) 



1 9 1 , r- 9® 

-mV 2 + V tI -- pA^+Gn + m— = 6 

2 2m \/n at 

where V is the condensate velocity and $ is the velocity potential, V = V$. The trapping potential plays the role of 
gravitational potential in hydrodynamics. We look for small perturbation from the equilibrium state. The equation 
for perturbation in the velocity potential for each of the component is the same as the one component case ||l8| 

^$ _ 2r Q 4- + (R 2 - r 2 )A$ = 0. (7) 
u>q or 

Continuity of the pressure at the interphase results in the following boundary condition 

LolRv^P- - LU 2 ®! =KU 2 R*^ "KW 2 $ 2 . (8) 

or Or 
Continuity of V r gives another boundary condition at the interface 

dr dr 

Eq. (0) and the boundary conditions (||), ^ constitute a complete set of equations necessary to determine the 
normal modes of the system. We expand the velocity potential in terms of spherical harmonics $ = $>(r)Yi m (0, <f>). 
From (^) the radial component can be written in terms of hypergeometric functions $1 = C\r l F{a, 0, 1 + 3/2, r 2 / Rf), 
= C 2 r l F(a, (3, 1, 1 — r 2 / R\) : where a + (3 — I + 3/2, a/3 = (I — uj 2 /uiq)/2. From matching boundary conditions, we 
finally obtain the eigenvalue equation for the normal mode frequencies 



to 2 [Z(7 + 3/2) + {l^uj 2 /aj 2 )Xx Sl (uj,x)] [l - {I ~ u 2 /u%)xs 2 (lo,x)] 

u 2 ~ (1 ~ K ' [1(1 + 3/2)(l - k) - xjl - uj 2 /lu 2 ) {k\ Si (uj, x) + {I + 3/2)s 2 (w, x))} ' 



where si (w, x) = F(a + 1,(3+1,1 + 5/2, \x)/F{a, (3, 1 + 3/2, Ax), s 2 (u, x) = F(a + 1,(3+1,2,1- x)/F(a, (3,1,1- x), 



x = Rl/Rl 
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One special solutions of Eq. ([lO|) is ui 2 = Iuj 2 , , which coincide with those for the one component condensate [Gl| . 
For this solution the components oscillate in-phase and $1 = k$2 oc r l Yi m (6, <fi). Another special exact solution is 
5u>Q with $1 a 1 - 5r 2 /3i? 2 , $2 oc A(l — 5r 2 /3R 2 ). For this solution the components oscillate out-of-phase if 



^Ji/hRi < R* < -\/3/5i?2 and in-phase otherwise. In general we have solved Eq. (10) numerically and found the 
normal mode frequencies u of the two component condensate as a function of the parameter k = \J a\\j a 22 for a 
fixed ratio R*/R 2 - The ratio -R*/i? 2 can be directly measured experimentally. Fig. 1 shows the low frequency modes 
that become imaginary at k > 1. These modes are peculiar for the two component systems and are analogous to the 
waves at the interphase between two layers of immiscible fluids under gravity jl9) . For this system, when the top layer 
becomes denser, the gravitational energy becomes higher and the system becomes unstable. In a similar manner as 
soon as k becomes greater than 1, our system becomes unstable. 

In the region |1 — k\ <C lthe mode frequencies are small: \ui\ <C u>q. We obtain 



LO 2 ^U 2 (l- K )f(l, X ), (11) 



where 



l(l + 3/2 + x ai (0,x))(xs 2 {0,x)-l) 
J(,X) z[ Sl (0,x) + (Z + 3/2)s 2 (0,z)] 

Eq. ( |Tl| ) describes the behavior of the low frequency modes in the region close to the point of instability k = 1. In 
this region oj oc \/l — k and becomes imaginary when k > 1. We next address the correction of this equation due to 
the finite overlap of the wave function and the surface tension. 

In the TF limit the interface tension results in small corrections of the order of £/i?2 to the normal mode frequencies. 
However, for the low frequency modes in the vicinity k « 1 the effect of interface tension is substantial because the 
mode frequencies themselves are close to zero. The interface tension a modifies the boundary condition for the 
pressure so that the pressure difference at the interface is equal to the surface pressure Pi — P2 = u{\jr\ + 1/V 2 ), 
where n, r 2 are the principle radii of curvature. As a result, the ratio of densities at the interphase is 




G 22 n\R* 



The interface tension is given by 



^^i+^Mv^^-m (12) 



where & represents the single condensate coherence length & = h/y / 2miGuni, the pressure P as Gun 2 / 2 and the 
condensate densities n, are estimated near the interface. Using the modified boundary condition Eq. ( |l l| ) is changed 
to 

9 , A 3cr ail - l)(l + 2) \ , 

\ 2R Jr P ran\ui^R% J 

The interface tension shifts the frequencies of the lowest modes and changes the stability region. For I — the inner 
droplet moves as a whole without changing its shape. However, the displacement of the droplet into the less dense 
region decreases the interface energy which is proportional to n 3 / 2 . This is the origin for the contribution —ia/R^P 
to the mode frequency. As k increases the mode with I — 1 becomes imaginary first. This determines the system's 
stability limit. The two component condensate is locally unstable when 

^ > 1 - = 1 - ^(WV^- D 1/2 . (14) 
Ct22 R*P oil* 

We next turn our attention to the calculation of the global stability condition. In the asymmetric phase, we have 
component one on the right and component two on the left (Fig. 2). In the TF approximation, the density distribution 
is again given by Eqs. (||), (^). The position of the interphase boundary between them is now different. We determine 
the global stability condition by equating the energy of the symmetric and the asymmetric phases. The calculation is 
similar to our previous work (l6[ . In the limit N\ -C N 2 we find a simple expression for the line of global instability: 
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(ai2/V a n a 22 - i) 1 



(15) 



In general, the calculation can only be carried out numerically. In Fig. 3 we plot the phase diagram that shows 
different stability regions of two component condensates. In our estimates we take = 0.01. When a\2 < y^oHoii 
the homogeneous binary mixture is a stable state. Otherwise the two components are phase separated. In the phase 
separated region for small ratios of ai 1/0122 the state (1,2) with the first component inside and the second outside is 
the only stable configuration. As the ratio 011/022 increases this configuration becomes globally unstable when we 
cross the left (solid) curve. However, the system is locally stable since the configuration corresponds to a local minima 
of energy. As 011/022 is further increased, we cross the line of local stability (second solid line) and the (1,2) state 
becomes unstable towards transforming into a new stable asymmetric state. In between the two solid lines, the system 
may tunnel quantum-mcchanically from the (1, 2) state to the asymmetric phase. A possible scenario is via quantum 
nucleation. How this takes place in detail is beyond the scope of the present paper. Eventually the asymmetric phase 
becomes the (2, 1) state as k is further increased. If initially the system is prepared in the (2, 1) state then as the 
ratio 011/022 is decreased this configuration first becomes globally unstable when we cross the right dotted line and 
locally unstable when 011/022 is close enough to 1. 

In summary, we have elucidated in detail a new symmetric-asymmetric transition in mixtures of BEC's. To simplify 
the numerical details, we have presented our results assuming the trapping frequencies of the two components are the 
same. Experimentally, the transition can be observed by either changing the interaction strengths or the ratio of the 
trapping frequencies. 
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Fig. 1. Low frequency modes as a function of k = -\/ a ii/ a 22 for different I = 1,2,3. The position of the interface is 
R, — R2/2. Frequencies become imaginary at k > 1. 
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Fig. 2. Symmetric and Asymmetric phases of two component BECs 
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Fig. 3. Phase diagram of the two component condensate in coordinates displaying relative interaction strength between 
bosons. 
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